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摘要 :【 目 的】 了解 小 红 珠 绢 蝶 Parnassius nomion 线粒体 基因 组 的 特征 ,并 从 线粒体 基因 组 水 平 探 
究 蝶 类 高 级 阶 元 的 系统 发 育 关系 。【 方 法 】 采 用 PCR 扩 增 技术 及 Sequencher 4. 8 拼接 软件 获得 小 
红 珠 纲 蝶 线 粒 体 基 因 组 全 序列 。 参 考 鲜 翅 目 昆 贝 已 知 线粒体 基因 组 全 序列 并 使 用 MEGA6.0 软件 
对 小 红 珠 绢 蝶 线 粒 体 基因 组 中 各 基因 进行 定位 和 注释 。 采 用 tRNA Scan-SE 1.21 在 线 预测 小 红 珠 
绢 蝶 线 粒 体 基因 组 tRNA 基因 的 二 级 结构 。 基 于 线粒体 基因 组 13 个 蛋白 质 编码 基因 的 核 苦 酸 序 
列 重建 了 包含 凤 蝶 总 科 中 凤 蝶 科 ( Papilionidae) 、 绢 蝶 科 ( Parnassiidae ) 、 粉 蝶 科 ( Pieridae ) 、 眼 蝶 科 
( Satyridae ) , 3k 9E ( Nymphalidae ) 、 灰 蝶 科 (Lycaenidae) , S£ 9k 4} ( Danaidae ) 、 珍 蝶 科 ( Acraeidae) , 
"f SEE ( Libyheidae ) fe YL 4} ( Riodinidae)10 个 科 28 种 蝴蝶 的 系统 发 育 关 系 。 【结果 】 结 果 表明 ， 
小 红 珠 绢 蝶 线 粒 体 基因 组 全 序列 总 长 度 为 15 362 bp (GenBank 登录 号 ; MF496134) , &,4- 13 个 蛋白 
质 编码 基因 、22 个 tRNA 基因 2 个 TRNA 基因 和 1 个 A+T 富 含 区 。 小 红 珠 绢 蝶 线 粒 体 基 因 组 中 存在 
较 高 的 A +T 含 量 (79.6%)。 小 红 珠 绢 蝶 线 粒 体 基 因 组 13 个 蛋白 质 编码 基因 中 UUA 的 相对 同 义 
密码 子 使 用 频率 (RSCU) 最 高 (5.08), 而 AGG 和 CCG 相对 同 义 密 码 子 使 用 频率 (RSCU ) 均 较 低 
(0), 3-5 K% Sasakia charonda coreana 的 分 析 结 果 一 致 。 在 所 测 得 的 22 个 tRNA 基因 中 , 除 
tRNA C99 缺少 DHU AE 2b, AA tRNA 基因 均 能 形成 典型 的 三 叶 草 结构 ,这 与 鲜 翅 目 中 目前 已 得 
到 的 其 他 昆虫 线粒体 基因 组 中 tRNA 基因 的 二 级 结构 一 致 。 系 统 发 育 分 析 结 果 显 示 , 凤 蝶 总 科 内 
遇 蝶 科 与 灰 蝶 科 的 亲缘 关系 最 近 ; 粉 蝶 科 与 映 蝶 科 、 珍 蝶 科 、 眼 炬 科 、 斑 蝶 科 、 吃 蝶 科 、 蝎 蝶 科 和 灰 蝶 
科 的 系统 发 育 关 系 更 近 ; 绢 蝶 科 与 凤 蝶 科 锯 凤 蝶 亚 科 亲缘 关系 最 近 , 随 后 二 者 与 凤 蝶 亚 科 物 种 聚 为 
一 支 。 在 绢 蝶 科 中 ,小 红 珠 绢 蝶 与 依 帕 绢 蝶 Parnassius epaphus 的 亲缘 关系 最 近 。【 结论 】 本 研究 支 
持 绢 蝶 科 物种 归 为 绢 蝶 亚 科 , 绢 蝶 亚 科 、 锯 凤 蝶 亚 科 和 凤 蝶 亚 科 归 入 凤 蝶 科 , 且 绢 蝶 亚 科 与 锯 凤 蝶 
亚 科 为 姐妹 群 。 

关键 词 : RA; 凤 炬 总 科 ; 小 红 珠 绢 蝶 ; 线粒体 基因 组 ; 系统 发 育 
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Abstract: [ Aim] The objective of this study is to understand the characteristics of the complete 
mitochondrial genome of Parnassius nomion and to explore the higher-level phylogenetic relationships of 
butterflies based on mitochondrial genome sequences. [ Methods] The PCR amplification technique and 
Sequencher 4. 8 software were used to obtain the whole sequence of mitochondrial genome of P. nomion. 
Based on the known complete mitochondrial genomes of lepidopteran species, MEGAG. 0 software was 
used to locate and annotate each gene of P. nomion mitochondrial genome. tRNA Scan-SE 1. 21 was used 
to predict the secondary structure of the tRNA genes online. Based on the nucleotide sequences of 13 
protein-coding genes of mitochondrial genome, the phylogenetic relationships among 28 butterflies of 10 
families ( Papilionidae, Parnassiidae, Pieridae, Satyridae, Nymphalidae, Lycaenidae, Danaidae, 
Acraeidae, Libytheidae and Riodinidae) of Papilionoidea were reconstructed. [Results] The results 
indicated that the complete mitochondrial genome of P. nomion is a circular molecule of 15 362 bp 
(GenBank accession no. : MF496134) , including 13 protein-coding genes, 22 tRNA genes, 2 rRNA 
genes and an A + T-rich region. The mitochondrial genome of P. nomion is biased toward a high A + T 
content (79. 6% ). 
protein-coding genes in the mitochondrial genome of P. nomion, the relative synonymous codon usage 
value is 5.08, but AGG and CCG has a relatively low RSCU of 0, which is consistent with the analysis 


results of Sasakia charonda coreana. All 22 tRNA genes show the classic clover-leaf structure, except 


The relative synonymous codon usage ( RSCU) of UUA is the highest in the 13 


that the dihydrouridine ( DHU) arm of tRNAS” ‘CY forms a simple loop, which is consistent with the 
secondary structure of tRNA genes in the mitochondrial genomes of other lepidopteran insects. The results 
of phylogenetic analysis showed that Riodinidae and Lycaenidae are sister groups, Pieridae is closer with 
the clade of ( Lycaenidae + Riodinidae) + ( ( ( Nymphalidae  Acraeidae) + Satyridae + Danaidae ) + 
Libyheidae) , Parnassiidae is a sister group of Zerynthiinae, and then Parnassiidae + Zerynthiinae is 
clustered with Papilioninae. P. nomion is closely related with Parnassius epaphus in Parnassiidae. 
[Conclusion ] Our study supports that Parnassiidae species should be included in Parnasiinae, 
Parnasiinae,, Zerynthiinae and Papilioninae all belong to Papilionidae, and Parnasiinae and Zerynthiinae 


are sister groups. 


Key words: Lepidoptera; Papilionoidea; Parnassius nomion; mitochondrial genome; phylogeny 
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线粒体 基因 组 ( mitochondrial genome, mtDNA 
genome ) 是 真 核 细胞 中 除 核 基 因 之 外 的 细胞 质 DNA 
遗传 信息 ,存在 于 线粒体 中 。 除 少数 较为 低 等 的 真 
核 生物 外 ,多 数 真 核 生 物 线粒体 基因 组 为 一 个 双 链 
闭合 的 环 状 DNA 分 子 ( 陈 念 和 赖 小 平 , 2011). 2X 
粒 体 基因 组 具有 结构 简单 稳定 、 相 对 分 子 量 较 小 、 核 
荫 酸 序列 比较 保守 不 含 或 仅 含 很 短 的 间隔 序列 和 
重复 序列 . 子 代 的 遗传 现象 与 母 本 性 状 一 致 以 及 点 
突变 率 高 且 便 于 在 较 短 的 时 间 内 检测 到 基因 变化 等 
特点 。 因 此 ,线粒体 DNA 被 广泛 用 于 生物 类 群 不 同 
分 类 水 平 的 系统 发 育 关系 、 种 群 遗传 结构 以 及 系统 
地 理学 等 方面 的 研究 (李靖 等 ,2010; ARES, 
2012; PSE, 2013) 。 近 年 来 ,利用 线粒体 基因 组 
全 序列 检测 和 分 析 物 种 间 分 子 系统 进化 的 方法 已 得 
到 广泛 关注 ,并 且 获 得 了 较为 重要 的 分 子 标记 信息 
( 王 菊 平等 ,2015; Sun et al., 2015; Hua et al., 





























2016; Liu and Zhao, 2016; Zhang et al., 2016; Chu 
et al., 2016; Z6, 2017) 。 

目前 ,在 国内 较为 通用 的 分 类 系统 中 ,小 红 珠 绢 
WE Parnassius nomion >Ë Jgi T 83H H ( Lepidoptera ) 绢 
TEE ( Parnassiidae ) 绢 蝶 属 Parnassius。 诸 多 国内 外 
专家 对 于 小 红 珠 绢 蝶 所 属 的 绢 蝶 科 的 分 类 地 位 仍 存 
在 较 大 分 歧 , 李 传 隆 . 宋 宝 云 和 周 阁 将 绢 蝶 、 风 蝶 作 
为 科 级 阶 元 处 理 ( 李 传 隆 和 朱 宝 云 ，1992; AF, 
1998) 。 李 传 隆 和 宋 宝 云 (1992 ) 认为 绢 蝶 科 包含 绢 
蝶 亚 科 Parnassiinae 和 亚 凤 蝶 亚 科 Zerynthiinae ( 包 
含 尾 凤 蝶 属 Bhutanitis, EHR JE Luehdorfia 和 丝 凤 
ME) Sericinus)  Jš|>ë ( 1998 ) 认为 绢 蝶 类 物种 应 为 
独立 的 科 级 阶 元 , 凤 蝶 科 包 含 凤 蝶 亚 科 Papilioninae 
和 锯 凤 蝶 亚 科 Zerynthiinae。 武 春生 将 绢 蝶 、 凤 蝶 、 
锯 凤 蝶 作 为 亚 科 级 阶 元 处 理 , 统 归 于 风 蝶 科 
Papilionidae, Jf. H. iA 2y 2H W& M. EE E A 28 9E Jm 
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Parnassius FIZ 2E Bë Jg Hypermnestra. Vj JB , 4 KUR VV. 
科 包 含 尾 风 蝶 属 Bhutanitis, Pë RUB 8 Luehdorfia 和 
ZZ ARJ Sericinus 3 属 ( 武 春生 ,，2001) 。 丁 亮 等 
(2007) 对 锯 凤 蝶 类 和 绢 蝶 类 的 分 类 地 位 及 系谱 关 
系 进行 了 研究 ,认为 绢 蝶 族 很 难 独 立 于 凤 蝶 科 而 独 
自 形成 科 级 分 类 单元 ,作为 亚 科 级 处 理 最 佳 ,并 且 认 
为 锯 凤 蝶 族 和 绢 蝶 族 亲缘 关系 十 分 密切 。 部 分 国外 
专家 将 绢 蝶 族 Parnassiini 与 锯 凤 蝶 族 Zerynthiini 作 
为 姐妹 群 , 同 为 亚 科 级 (Munroe，1961; Hancock, 
1983; Igarashi, 1984) 。 

鉴于 近年 越 来 越 多 的 研究 人 员 采 用 线粒体 多 基 
因 联 合 序列 构建 物种 的 系统 发 育 树 来 探究 各 物种 之 
间 的 系统 发 育 关 系 和 分 类 地 位 ( 杨 孝 文 ,， 2015; 郭 
EAE, 2017; HE, 2017) ,本 研究 测定 了 小 红 
珠 绢 蝶 线 粒 体 基因 组 全 序列 并 在 GenBank 中 下 载 
获得 与 小 红 珠 绢 蝶 相 关 的 28 种 包含 绢 蝶 类 、 凤 蝶 
类 、 粉 蝶 类 等 物种 的 线粒体 全 基因 组 序列 ,采用 最 大 
简约 法 构建 其 系统 发 育 树 ,探究 小 红 珠 绢 蝶 所 在 的 
高 级 阶 元 的 分 类 地 位 ,以 期 为 绢 蝶 科 物种 的 分 类 地 
位 研究 提供 新 的 基础 知识 ,并 溢 清 绢 蝶 类 各 物种 之 
间 的 系统 发 育 学 关系 。 








































































































1 材料 与 方法 


1.1 实验 材料 

本 研究 所 用 的 小 红 珠 绢 蝶 个 体 采 自 中 国 山西 省 
五 台 县 (海拔 2 795 m, 平 均 气 温 20%C ) ,野外 所 采集 
的 小 红 珠 绢 蝶 成 虫 活体 标本 置 于 三 角 纸 袋 中 保存 ， 
并 于 24 h 内 冻 存 于 -800 冰箱 中 备用 。 
1.2 基因 组 DNA 的 提取 和 纯化 

将 小 红 珠 绢 蝶 标 本 从 -80% 冰箱 中 取出 后 置 于 
冰 上 ,选取 发 育 良好 的 个 体 分 离 其 胸部 肌肉 ,随后 采 
用 OMEGA Insect DNA Kit 试剂 盒 并 参考 酚 / 氧 仿 抽 
提 样 本 总 DNA 的 方法 进行 小 红 珠 绢 蝶 样 本 总 DNA 
的 提取 。 采 用 1. 0% 琼脂 糖 凝 胶 电 泳 检测 其 浓度 ， 
使 用 分 光 光 度 计 推 测 其 浓度 ,后 置 于 -20%C 保存 
备用 。 
1.3 线粒体 基因 组 全 序列 的 PCR 扩 增 和 测序 

采用 Clustal X 1. 8 软件 (Thompson et al., 
1997) 对 GenBank 下 载 的 部 分 相关 昆虫 线粒体 基因 
组 全 序列 进行 比 对 , 找 出 同 源 性 较 高 的 保守 序列 。 
针对 这 些 保守 区 序列 使 用 Primer Premier 5.0 软件 
(Singh et al., 1998 ) 设计 多 对 涵盖 小 红 珠 绢 蝶 线 粒 
体 基因 组 全 序列 的 PCR 引物 ,对 于 采用 上 述 引 物 未 












































扩 增 出 的 序列 ,我 们 结合 已 经 获得 的 序列 采用 
Primer Premier 5.0 软件 重新 设计 小 红 珠 绢 蝶 的 特异 
引物 进行 PCR 扩 增 ( 表 1)。 使 用 设计 好 的 引物 和 
TaKaRa LA Tag 聚合 酶 进行 小 红 珠 绢 蝶 线 粒 体 基因 
组 序列 片段 的 PCR 扩 增 。 根据 TaKaRa LA Taq 5€ 
合 酶 提供 的 使 用 说 明 设计 PCR 扩 增 体系 (25 uL): 
DNA 模板 (25 ng/uL)1 uL, 10 x LA PCR Buffer 2.5 
L， 上 下 游 引物 (10 umol/L) 44 1 uL, dNTPs (2.5 
mmol/L)4 uL, TaKaRa LA Taq 酶 (5 U/uL) 0.25 
pL, X 8 zk P E $8] 25 uL, PCR 扩 增 反应 条 件 : 
95*C 预 变性 5 min; 94% 变性 30 s, 50 -60% 退火 30 
s, 72% 延伸 1 min, 共 35 个 循环 后 ;再 72% 延伸 10 
min, PCR 产物 使 用 1% 琼脂 糖 凝 胺 电泳 检测 。 用 
OMEGA 公司 (美国 ) 胶 回收 试剂 盒 回 收 和 纯化 目的 
条 带 后 转 人 大 肠 杆 菌 Escherichia coli DH 5a 感受 态 
细胞 ,筛选 阳性 克隆 并 检测 ,委托 上 海 立 菲 生物 技术 
有 限 公 司 进行 双向 测序 。 
1.4 全 序列 的 拼接 .分 析 与 注释 

对 得 到 的 DNA 序列 片段 采用 Clustal X 1.8 软 
件 与 设计 引物 所 用 序列 进行 比 对 ,将 所 有 比 对 正确 
的 序列 置 人 软件 Sequencher 4. 8 ( Gene Codes 
Corporation, 2008) 中 进行 拼接 ,得 到 小 红 珠 绢 蝶 线 
粒 体 基 因 组 全 序列 ,并 使 用 在 线 序列 提交 软件 
Bankit 将 其 全 序列 提交 GenBank 中 , 登录 号 为 
MF496134。 参 考 鳞 翅 目 昆虫 中 已 知 线粒体 基因 组 
全 序列 对 小 红 珠 绢 蝶 的 线粒体 基因 组 中 各 基因 进行 
定位 和 注释 ( Kim SR er al., 2009; Jiang et al., 
2009) ,使 用 ClustalX 2. 1 软件 对 小 红 珠 绢 蝶 和 相关 
鳞 翅 目 昆 虫 线粒体 基因 组 全 序列 进行 比 对 ,确定 基 
因 的 位 置 。 使 用 MEGA 6 软件 计算 小 红 珠 绢 屿 线 粒 
体 全 序列 的 核 昔 酸 组 成 ,包括 A+T 含 量 ,C + E 
量 等 ,并 分 析 和 蛋白 质 编码 基因 的 氨基 酸 使 用 情况 及 
相对 同 义 密 码 子 使 用 频率 (relative synonymous 
codon usage, RSCU), 。 采 用 tRNA Scan-SE 1. 21 
(http: //lowelab. ucsc. edu/tRNA Scan-SE ) 在 线 预 
测 小 红 珠 绢 蝶 线 粒 体 基因 组 tRNA 基因 的 二 级 结构 。 
1.5 系统 发 育 分 析 

基于 线粒体 基因 组 和 蛋白质 编码 基因 的 核 苷 酸 序 
列 ,采用 最 大 简约 法 (Maximum Parsimony method, 
MP) 和 贝 叶 斯 法 (Bayesian inference, BI) 以 美国 白 
蛾 Hyphantria cunea, Zum Bombyx mori 和 亚洲 玉米 
Bs Ostrinia furnacalis 作为 外 群 ,构建 了 弄 蝶 总 科 和 
JA ES ËL pg QÇ P& BE (Papilionidae ) , 28 W 科 
( Parnassiidae ) 、 粉 蝶 科 (Pieridae ) IRIRE} Satyridae) , 
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表 1 本 研究 使 用 的 PCR 引物 
Table 1 Primers used in this study 




















引物 引物 序列 (5” -37) 片段 大 小 (bp) 参考 文献 
Primers Nucleotide sequence Product size References 
TM-J210 AATTAAGCTATTAGGTTCATACCC Simon et al., 2006 
TY-N1433 GGCTGAATAATAAGCGATAAATTGTAAA 1223 Simon et al., 2006 
TW-J1301 GTTAAAWTAAACTAATARTCTTCAAA Simon et al., 2006 
C1-N2353 GCTCGTGTATCAATATCTATWCC 1 052 Simon et al., 2006 
C1-J2195 TGATTTTTTGGWCATCCWGAAGT Simon et al., 2006 
C2-N2776 GATAATCTGAATATCGWCGNGG 581 Simon et al., 2006 
C1-J2756 ACATTTTTTCCACARCAYTT Simon et al., 2006 
C2-N3389 TATTCATARCTTCARTATCATTG 633 Simon et al., 2006 
TL2-J3043 GGCAGACTATATGYAATGRATTTAA 本 研究 This study 
TK-N3796 ACTATAAAATGGTTTAGAG 753 本 研究 This study 
A6-]4463 TTTATTCATATWATWCCNCAAGG Simon et al., 2006 
C3-N4908 CGTGAAAYATCTCGTCATCATTG 763 Simon et al., 2006 
C3-N4725 TCGACCTGGAACTTTAGC Hong et al., 2009 
TN-N6354 'TGGATCAAATCCACATTCA 445 Hong et al., 2009 
TN-J6172 AGAGGTATATCACTGTTAATGA Simon et al., 2006 
N5-N7211 TTAAAGCTTTAYTATTTATRTGYGC 939 Simon et al., 2006 
N5-]7077 TTAAATCTTTWGARTAAAAYCC Simon et al., 2006 
N5-N7793 TTGGGTTGRGATGGNYTAGG 914 Simon et al., 2006 
N5-]7572 AAAAGGAATTTGAGCACTTTTWGT Simon et al., 2006 
N4-N8727 AAATCTTTRATTGCTTATTCWTC 1 039 Simon et al., 2006 
N4-N7657 CACCAGCAGTAACTAATGTTGATG Simon et al., 2006 
N4-N7967 CAACCTGAACGAATACAAGCTGG 716 Simon et al., 2006 
NA-N7817 GGAGCTTCTACATGAGCTTTTGG Hong et al., 2009 
N6-N9967 GACCTCAAGGTAAAACATAACC 1155 Hong et al., 2009 
N6-N9867 TATCACAACCTAACAT 本 研究 This study 
N6-N10067 AATGAGATATTAGATGGGGAAGG 988 本 研究 This study 
N6-N9954 CGAACTATTCACGCTAATGGAGC Hong et al., 2009 
N6-N11967 CTTATCGAYAAAAAAGWTTGCGACCTCGATGTTG 1 838 Hong et al., 2009 
N1-711876 CGAGGTAAAGTMCCWCGAACYCA Simon et al., 2006 
N1-N12595 GTWGCTTTTTTGACTTTATTRGARCG 1 446 Simon et al., 2006 
N1 -J12261 TACTTCATAAGAAATAGTYTGRGC Simon et al., 2006 
LR-N13000 TTACCTTAGGGATAACAGCGTAA 719 Simon et al., 2006 
LR-J12888 CCGGTTTGAACTCARATCATGTA Simon et al., 2006 
LR-N13889 ATTTATTGTATCTTKTGTATCAG 739 Simon et al., 2006 
LR-N13779 CGGTTTGAACTCAGATCATGTAAG Hong et al., 2009 
LR-N14089 TATTGTATCTTGTGTATCAGAGTTTA 1 101 Hong et al., 2009 
LR-J13342 CCTTTGTACRGTCAAAATACYGC Simon et al., 2006 
SR-N14220 ATATGYACAYATTGCCCGTC 878 Simon et al., 2006 
SR-J14197 TACCYCTACTTTGTTACGACTT Simon et al., 2006 
SR-N14745 GTGCCAGCAGYYGCGGTTANAC 548 Simon et al., 2006 
SR-J14610 ATAATAGGGTATCTAATCCTAGT 330 Simon et al., 2006 
TM-N200 TCCTTTATATRTGAGGTATGARCC Simon et al., 2006 
SR-J14750 TGCTGGCACAAAATTAGTTA 358 本 研究 This study 
TM-N280 AATCCAATTCAGCAACCTAA 本 研究 This study 
FW2018 TTCTATGAAGATTAGGATTTGTTTT T 本 研究 This study 
RV3030 AATGAATTACGTCAGTAGCAGTTAC 本 研究 This study 
FW8108 AAACCAAAAGAAATTAAATCACAAC dé 本 研究 This study 
RV9104 GATAGACGTAAAGGACCAAAAAAAA 本 研究 This study 
FW10152 TGAAGACCCTTATATTATTACTGGT Sia 本 研究 This study 
RV11126 CTTTTCTTACATTATTAGAACGGAA 本 研究 This study 
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WIRE RE Nymphalidae) 、 灰 蝶 科 ( Lycaenidae ) REIRE} 
( Danaidae ). 32 W ËL ( Acraeidae ), WA HR Æ} 
( Libyheidae) # BL £|. ( Riodinidae) 10 A (J8]38, 








1998) f 28 种 蝴蝶 (包括 小 红 珠 绢 蝶 ) 的 系统 发 育 MAR. 


表 2 系统 发 育 分 析 所 用 类 群 信息 


行 序列 比 对 后 ,将 序列 上 








列表 


Table 2 List of taxa used for phylogenetic analysis 





树 ( 表 2)。 将 上 述 28 种 蝴蝶 及 外 群 的 线粒体 蛋白 
编码 基因 核 苷 酸 序列 分 别 采 用 ClustalX 2. 1 软件 进 
EKEK, 以 防 出 现 移 码 


GenBank 登录 号 













































































































































































总 科 名 科 名 亚 科 名 物种 
Superfamily Family Subfamily Species GenBank 
accession number 
外 群 TENURE Bombycidae 2 Bombyx mori NC. 002355.1 
Outgroup IRIRE} Pyralidae NPM CM Ostrinia furnacalis NC_003368 
灯 蛾 科 Arctiidae 美国 白 峨 Hyphantria cunea NC_014058 
弄 蝶 总 科 JEME} Hesperiidea FEWE IE} Hesperiinae 卵 链 弄 蝶 Heteropterus morpheus NC_028506 
Hesperoidea FAJI Parnara guttata NC_029136 
凤 蝶 总 科 粉 蝶 科 Pieridae 黄粉 蝶 亚 科 Coliadinae 宽 边 黄粉 蝶 Eurema hecabe KF. 881047 
Papilionoidea 粉 蝶 亚 科 Pierinae 黑 纹 粉 蝶 Pieris melete NC_010568 
袖 粉 蝶 亚 科 Dismorphiinae 莫 氏 小 粉 蝶 Leptidea morsei NC_022686 
眼 蝶 科 Satyridae 锯 眼 蝶 亚 科 Elymninae HM E Lethe dura KF_906485 
腿 蝶 亚 科 Satyrinae All s Callerebia suroia NC. 026060 
亚洲 白眼 蝶 Melanargia asiatica KF_906486 
峡 蝶 科 Nymphalidae Digit Ip |. Apaturinae Me DJ BEBE Apatura ilia JF_437925 
JK BERE. Lycaenidae £g JK BE VV E|: Theclinae 黄 灰 蝶 Japonica lutea NC_026517 
眼 灰 蝶 亚 科 Polyommatinae Cupido argiades NC_020779 
灰 蝶 亚 科 Lycaeninae 红 灰 蝶 Lycaena phlaeas NC_023087 
绢 蝶 科 Parnassiidae 依 帕 绢 蝶 Parnassius epaphus NC_026864 
B EHHE Parnassius imperator NC_025587 
元 首 绢 蝶 Parnassius cephalus NC_026457 
红 珠 绢 蝶 Parnassius bremeri NC_014053 
fof iic 89 28 8E Parnassius apollo NC. 024727 
小 红 珠 绢 蝶 Parnassius nomion This sthdy 
JAEP} Papilionidae 锯 凤 蝶 亚 科 Zerynthiinae AUR UE Luehdorfia taibai NC_023938 
丝带 凤 蝶 Sericinus montela HQ_259122 
凤 蝶 亚 科 Papilioninae Papilio glaucus NC. 027252 
JE XUME Papilio protenor NC. 034317 
宽 尾 凤 蝶 Papilio maraho NC_014055 
A BEI UM Leptidea morsei NC, 014398 
珍 蝶 科 Acraeidae +R Acraea issoria NC. 013604 
WEBER Libyheidae F IRI Libythea. celtis NC. 016724 
JL EL Riodinidae HIRIE Abisara fylloides NC. 021746 
斑 蝶 科 Danaidae ABE Danaus chrysippus KF. 690637 


MP 分 析 在 PAUP ° 4b4a 软件 中 完成 ,选择 启发 
式 搜 索 (heuristic search ) 、 树 二 等 分 再 连接 (tree- Uu 
bisection reconnection, TBR) 和 100 次 随机 加 入 序列 
(addition sequence replicates) 等 参数 来 搜索 MP 树 ， 
同时 采用 了 非 参 数 自 展 法 (non-parametric bootstrap ) 











1 000 次 重复 检验 各 节点 的 置信 和 度 。 

















分 析 采 用 MrBayes 3. 1. 2( Huelsenbeck , 


2001) 软件 来 完成 ,系统 发 育 树 中 各 分 支 的 可 信和 度 
用 后 验 概率 来 表示 (posterior probability, PP), XH 
Modeltest 3. 7 软件 (Posada and Crandall, 1998) 中 的 
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等 级 似 然 率 (hLRT) 和 AIC 检验 标准 估算 出 最 适合 
本 研究 中 数据 集 的 蔡 代 模型 , 即 "CTR +1+G”。 随 
后 在 MrBayes 3. 1.2 软件 中 ,采用 马尔 科 夫 蒙特 卡 
罗 方 法 (Markov Chain Monte Carlo, MCMC) fA Œ bl 
叶 斯 树 。 




















2 结果 


2.1 小 红 珠 绢 蝶 线 粒 体 基因 组 结构 

小 红 珠 绢 蝶 线 粒 体 基因 组 全 序列 总 长 度 为 
15 362 bp, 其 基因 组 成 与 其 他 已 知 鳞 翅 目 昆 虫 线 粒 
体 基因 组 成 完全 一 致 ,都 包含 了 后 生动 物 线粒体 基 
因 组 成 中 的 37 种 基因 , 即 13 个 蛋白 质 编码 基因 、22 
个 tRNA 基因 、2 个 rRNA XE ES ( 16S/rraL 和 128/ 
rs) 和 1 个 A+T 富 含 区 (A+T-rich region) (Liu et 
al., 2008; Yang et al., 2009) 。37 个 基因 中 的 9 个 
和 蛋白质 编码 基因 和 14 个 tRNA 基因 位 于 工 链 ;其 余 
4 个 蛋白 质 编码 基因 和 8 个 tRNA JE DSL f J H 
链 , 且 上 链 和 日 链 上 基因 编码 方向 相反 (图 1)。 小 
红 珠 绢 蝶 线 粒 体 基因 组 中 基因 重奏 区 共 13 Ab EA 
序列 最 长 的 区 域 位 于 tRNA 和 tRNA™ 基 因 之 间 ， 
共 8 bp; 基 因 间 隔 区 共 12 处 ,其 中 间隔 序列 最 长 区 
域 位 于 tRNA 和 ND2 基因 之 间 , 长 度 为 40 bp; BE 
无 重合 又 无 间隔 的 基因 11 处 ( 表 3)。 
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小 红 珠 绢 蝶 
Parnassius nomion 
15 362 bp 

















图 1 小 红 珠 绢 蝶 线 粒 体 基因 组 结构 
Fig. 1 





Structure of the mitochondrial genome 


of Parnassius nomion 


2.2 小 红 珠 绢 蝶 线 粒 体 基 因 组 核 背 酸 组 成 
采用 MEGA 6.0 软件 分 别 对 小 红 珠 绢 蝶 线 粒 体 
基因 组 的 全 序列 、13 个 蛋白 质 编码 基因 序列 .2 个 











rRNA 基因 和 A +T 富 含 区 序列 的 长 度 和 A+T 碱 基 
百 分 含 量 进行 统计 分 析 ( 表 4) ,结果 显示 ,小 红 珠 绢 
蝶 线 粒 体 基因 组 全 序列 中 A, T, G 和 C 含量 分 别 为 
39. 6% , 40% , 8. 296 和 12.2%,A+T 含 量 达 到 
79. 6% ,G + C 含量 为 20.4% ;13 个 蛋白 质 编码 基因 
序列 中 A + 了 含量 达到 80. 2%,G + C 会 量 为 
19. 8% ;rrnL 基因 序列 中 A + T 4 EU 84.496 , G + 
C 含量 为 15. 6% ; rra$ 基因 序列 中 A + T 8E? 
85. 2% ,G +C 含量 为 14.8% ;A+T 富 含 区 中 A +T 
含量 则 高 达 93% ,G + C 含量 为 7% 。 此 外 ,小 红 珠 
绢 蝶 线 粒 体 基因 组 全 序列 的 AT Aii SURE GC 偏 斜 分 
别 为 -0.005 和 -0. 196 ,这 表明 其 序列 中 有 较 多 的 
T fli C, 

将 小 红 珠 绢 蝶 碱 基 组 成 分 析 结 果 与 已 报道 的 
13 种 鳞 翅 目 昆 虫 线粒体 基因 组 全 序列 的 碱 基 含 量 
进行 比较 分 析 ( 表 5) ,结果 显示 ,小 红 珠 绢 蝶 线 粒 体 
基因 组 序列 的 A + T & = (79. 6% ) 略 高 于 仁 眼 蝶 
(79.1% ) ,但 比美 国 白 蛾 等 另外 12 PERA H Fih 
为 低 。 除 朝 灰 蝶 线 粒 体 基因 组 13 个 蛋白 质 编码 基 
因 序 列 的 A+T 含 量 较 高 (81.4% ) 外 ,小 红 珠 绢 蝶 
的 A+T 含 量 (80.2% ) 高 于 GenBank 数据 库 中 下 载 
的 其 余 12 种 鳞 翅 目 昆 虫 的 A+T 了 含量。 此 外 ,小 红 
珠 绢 蝶 A+ 了 富 含 区 的 A +T RER 93.0% ,高 于 
RE %* Antheraea pernyi (89. 596) , KÆ Antheraea 
yamamai(90. 4% ) 、 黑 纹 粉 蝶 Pieris melete( 89.296 ) 、 
KA edit (01. 696 ) ( Liu et al., 2008; Kim SR et al., 
2009; Hong et al., 2009; 王 菊 平 等 , 2012) 。 该 结果 
符合 鳞 翅 目 昆虫 线粒体 基因 组 具有 A + T (im ln] PEB 
寺 征 。 
2.3 蛋白质 编码 基因 
2.3.1 小 红 珠 绢 蝶 线 粒 体 基因 组 蛋白 质 编码 基因 
分 析 :13 个 蛋白 质 编码 基因 分 别 为 3 个 细胞 色素 C 
还 原 酶 亚 单位 复合 体 编码 基因 (C01, COIT 和 
COIII) ,ATP 酶 的 2 个 亚 基 基因 (47P6 和 47P8 ) ， 
NADH 脱氧 酶 的 7 个 亚 基 基因 (ND1，NVD2 ，VD3 , 
ND4, ND4L, ND5 和 ND6) 及 细胞 色素 b 基因 
(CyiB)。 这 13 个 蛋白 质 编码 基因 与 已 测 出 的 其 他 
鳞 翅 目 昆 虫 线粒体 基因 组 蛋白 质 编 码 基 因 的 排列 顺 
序 完 全 一 致 (Yukuhiro et al., 2002; Lee et al., 2006; 
Cameron et al., 2008) , 
2.3.2. 蛋白质 编码 基因 密码 子 使 用 情况 分 析 :13 
个 蛋白 质 编 码 基 因 中 , ND2, ATP8, ND3, NDS 和 
ND6 这 5 个 基因 以 ATT 作为 起 始 密 码 子 ; COIT, 
ATP6, COII 和 CytB 以 ATG 作为 起 始 密码 子 ; NDA 
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和 DN4L VA TTA 作为 起 始 密码 子 ;而 CO 以 CGA 作 
为 起 始 密码 子 ( 表 3 ) ,这 与 鳞 翅 目 中 其 他 物种 ,例如 
KA EE IESENE Areyreus hyperbius | Sl #8 V) I die 
Apatura ilia ( Wang et al., 2011; E3 Æ, 2012; 
Chen et al., 2012) 等 昆虫 线粒体 COI 基因 以 CGA 
为 起 始 密码 子 的 现象 一 致 。 此 外 , 小 红 珠 绢 蝶 线 粒 
体 除 ND3 基因 外 ,VD2, ATP8, ATP6, COII, NDS, 

















ND6 和 CytB 基因 均 以 标准 密码 子 TAA. 作为 终止 密 
码 子 的 现象 与 鳞 怒 目 其 他 物种 一 致 ,例如 大 紫 峡 蝶 
ND2, ATP8, ATP6, COIII, ND3, ND5, ND6 和 CytB 
基因 以 标准 密码 子 TAA 作为 终止 密码 子 ( 王 菊 平 
等 , 2012) ;而 COI, COH, NDA, ND4L 和 ND1 基因 
均 以 单个 T 作为 终止 密码 子 。 





RI 小 红 珠 绢 蝶 线粒体 基因 组 注释 结果 


Table3 Annotation of the mitochondrial genome of Parnassius nomion 


























基因 方向 位 置 (bp) 大 小 (pp) 反 密 码 子 起 始 密码 子 。 终止 密码 子 基因 间 间 隔 (bp) 
Gene Direction Location Size Anticodon Start codon Stop codon Intergenic nucleotides 
tRNAMet F 1 -66 66 CAT 
tRNA! F 67 -130 64 TAT 0 
tRNA" R 128 - 196 69 TTG -3 
ND2 F 231 -1 250 1 014 ATT TAA 40 
tRNA T? F 1 250 - 1315 66 TCA -1 
tRNASY R 1 308 -1 373 66 GCA -8 
tRNA T: R 1378-1441 64 GTA 4 
COI F 1 444 -2 975 1 532 CGA T 2 
tRNA Le (Uf) F 2975 -3 041 67 TAA -1 
COII F 3 042 -3 723 682 ATG T 0 
tRNA! F 3 724 -3 794 71 CTT 0 
tRNA ^» F 3 794 -3 860 67 GTC -1 
47P8 F 3 861 -4 025 165 ATT TAA 0 
ATP6 F 4 019 -4 696 678 ATG TAA -7 
COIII F 4 696 —5 484 789 ATG TAA -1 
tRNA GY F 5 488 -5 554 67 TCC 3 
ND3 F 5 555 -5 908 354 ATT TAA 0 
tRNA AN F 5 908 -5 973 66 TGC -1 
tRNA F 5 973 -6 038 66 TCG -1 
tRNA ^^ F 6 039 -6 102 64 GTT 0 
tRNASe(AGN) F 6 106 -6 163 58 GCT 3 
tRNA CI F 6 185 -6 250 66 TTC 21 
tRNAP™ R 6 249 -6 314 66 GAA -2 
NDS R 6 317 -8 049 1 733 ATT TAA 2 
tRNA R 8 050 -8 113 64 GTG 0 
ND4 R 8 113 -9 453 1 341 TTA T -1 
ND4L R 9 453 -9 743 291 TTA T -1 
tRNA r F 9 749 -9 812 67 TGT 5 
tRNA?” R 9 813 -9 877 65 TGG 0 
ND6 F 9 880 — 10 410 531 ATT TAA 2 
CytB F 10 428 - 11 579 1 152 ATG TAA 17 
tRNASerCUCN) F 11 578 - 11 644 67 TGA E) 
NDI R 11 661 - 12 599 939 CTA T 16 
tRNALe"(CUN) R 12 601 - 12 669 69 TAG 1 
rm, R 12 670 - 14 017 1 348 0 
tRNA VI R 14 018 - 14 082 65 TAC 0 
rrnS R 14 083 - 14 879 791 0 
控制 区 Control region 14 880 - 15 362 483 0 
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表 4 小 红 珠 绢 蝶 线 粒 体 基因 组 核 背 酸 组 成 
Table 4 Nucleotide composition of the complete mitochondrial genome of Parnassius nomion 

基因 序列 A T G C A+T G +C AT 偏 斜 GC 偏 斜 

Gene sequences (96) (96) (96) (96) (96) (96) AT-skew GC-skew 

全 基因 组 Whole genome 39.6 40 8.2 12.2 79.6 20.4 —- 0. 005 - 0.196 

蛋白 质 编码 基因 Protein-coding genes 34.4 45.8 10.1 9.7 80.2 19.8 —0.142 0.02 
rRNA 基因 (rmL) rRNA genes ( rrnL) 41.7 42.7 5.0 10.6 84.4 15.6 -0.012 -0.359 
rRNA 基因 (rrmS) rRNA genes (rrnS) 41.8 43.4 5.0 9.8 85.2 14.8 -0.019 -0.324 
A+T HAE A + T-rich region 44.3 48.7 2:4 4.3 93.0 7.0 —- 0.047 -0.229 
X5 己 知 的 鳞 翅 目 昆 虫 物 种 线粒体 基因 组 的 核 音 酸 组 成 
Table 5 Nucleotide composition of the known mitochondrial genomes of lepidopteran species 
PCG rrnL rrnS AAT RAR 
A + T-rich region GenBank 
物种 大 小 A+T 登录 号 

Species Size (96) 大 小 A+T 大 小 A+T 大 小 A+T 大 小 AET GenBank 

(bp) Size (9) Size (96) Size (96) Size (ai. Arcsin ño: 
(bp) (bp) (bp) (bp) 

XE HIR Hyphantria cunea 15481 80.4 10851 78.5 1426 85.0 808 84.5 357 95.0 NC. 014058 
REE Phthonandria atrilineata 15499 81.0 10849 79.1 1400 85.8 803 86.3 457 98.2 NC 010522 
FEX Antheraea pernyi 15 566 80.1 10875 78.4 1380 84.0 7716 84.4 334 89.5 NC_004622 
KRÆ Antheraea yamamai 15338 80.3 10863 78.8 1369 83.8 TIS 84.1 552 90.4 NC 012739 

A C Manduca sexta 15516 81.8 10831 80.1 1 391 85.3 TTI 85.7 324 95.4 NC. 010266 
RUNA- US HB, Grapholita molesta 15716 80.9 10901 78.9 1 377 84.8 712 85.3 TIO 96.0 NC 014806 
RAR Acraea issoria 15245 79.7 11186 78.1 1 331 83.9 788 83.7 430 96.0 NC. 013604 
仁 眼 蝶 Hipparchia autonoe 15489 79.1 10838 76.9 1 335 83.7 TIS 85.2 678 94.6 NC 014587 
黑 纹 粉 蝶 Pieris melete 15140 79.8 10826 78.5 1319 83.4 TTI 85.5 351 89.2 NC. 010568 
HIK! Coreana raphaelis 15314 82.7 10797 81.4 1330 85.3 TTI 85.8 375 94.2 NC 007976 
红 珠 绢 蝶 Parnassius bremeri 15389 81.3 10845 80.1 1344 83.9 773 85.1 504 93.6 NC_014053 
金 斑 喉 风 蝶 Teinopalpus aureus 15242 79.8 10834 78.2 1320 82.4 781 85.6 395 93.1 NC 014398 
KA k Sasakia charonda coreana 15233 79.9 11231 78.4 ] 330 84.0 TI6 84.9 381 91.6 JX. 119051 
小 红 珠 绢 蝶 Parnassius nomion 15362 79.6 11199 80.2 1 348 84.4 797 85.2 483 93.0 This study 




















PCG: 蛋白 质 编码 提 


本 文 对 小 红 珠 绢 蝶 线 粒 体 基因 组 的 13 个 蛋白 
质 编码 基因 的 氨基 酸 使 用 情况 (图 2: A) 及 相对 同 
义 密码 子 使 用 频率 RSCU 进行 统计 分 析 ( 图 2: B), 
13 个 蛋白 质 编码 基因 中 使 用 较 频 繁 的 氨基 酸 为 
Leu, He, Phe, Met 和 Ser; 使 用 相对 较 少 的 氨基 酸 
73 His, Asp, Gln, Arg 和 Cys( 图 2: A)。13 个 蛋白 
质 编码 基因 使 用 频率 相对 较 高 的 密码 子 为 UUA, 
AUU 和 UUU ,其 中 UUA 的 相对 同 义 密码 子 使 用 频 
率 最 高 ,其 RSCU 值 为 5.08 ,而 ACC 和 CCG 相对 同 
义 密码 子 使 用 频率 均 较 低 , 其 RSCU 值 为 0, 这 与 大 
KERIT R ER FE, 2012 ) 一 致 。 
2.4 tRNA 基因 及 结构 特征 

通过 与 鳞 运 目 其 他 昆虫 线粒体 基因 组 进行 同 源 
性 比 对 并 参考 tRNAScan-SE 在 线 预 测 结果 , 获得 小 
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因 Protein coding gene. 







































































红 珠 绢 蝶 P. nomion 线粒体 基因 组 22 个 tRNA 基因 
序列 。 序 列 分 析 结 果 显 示 ( 表 3) ,这 22 个 tRNA 基 
因 序 列 总 长 度 为 1 450 bp ,其 排列 顺序 与 除 蝙蝠 峨 
科 以 外 的 鳞 翅 目 昆 虫 完全 相同 。 小 红 珠 绢 蝶 线 粒 体 
基因 组 中 每 个 tRNA 基因 的 长 度 范围 均 在 58 bp 
(tRNA ) ~ 71 bp (tRNA”) 之 间 。 其 中 
tRNA'"", tRNA™, tRNA™, iRNA'^UUP ，tRNAD ， 
tRNA^®, tRNA (RNA^^,. (RNA tRNA", 
tRNA (RNA " (RNA T 和 iRNA C™ 14 个 
tRNA 基因 由 工 链 编码 , 男 外 tRNA, tRNA, 
tRNA™Y, tRNA®™, iRNA"5, tRNA™ 和 tRNA 
和 tRNA""8 个 tRNA 基因 由 H 链 编 码 。 

基于 tRNA 结构 特征 和 tRNAScan-SE 的 预测 结 
， 预 测 出 小 红 珠 绢 蝶 线 粒 体 基因 组 22 个 tRNA 基 
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Leu(L) Ile(I) Phe(F) Met(M) Ser(S) Asn(N) Tyr(Y) Gly(G) 
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) Lys(K) Pro(P) Trp(W) Glu(E) His(H) Asp(D) Gln(Q) Arg(R) Cys(C 
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Leu(L) Ile(I Phe(F) Met(M) Ser(S) Asn(N) Tyr(Y) Gly(G) 


) Thr(T) Val( Mia Ala(A) Lys(K) Pro(P) Trp(W) Glu(E) His(H 


) Asp(D) GIn(Q) Arg(R) Cys(C) 














图 2 小 红 珠 绢 蝶 线 粒 体 基因 组 13 个 蛋白 质 编码 基因 氨基 酸 的 使 用 频率 (A) 和 同 义 密码 子 相对 使 用 度 (B) 


Fig. 2 Usage of amino acids ( A) and the relative synonymous codon usage (B) of 13 protein-coding genes of 











mitochondrial genome in Parnassius nomion 


因 的 二 级 结构 (图 3) ,转运 氨基 酸 Leucine 和 Serine 
的 tRNA 分 别 对 应 两 个 tRNA 基因 , 即 IRNA DY, 

tRNAILen(CUN) 和 tRNA SAN) : tRNASer(UCN) , 转运 运 其 余 
20 个 氨基 酸 的 tRNA 各 对 应 一 个 tRNA 基因 。 分 析 
这 22 个 tRNA 基因 的 二 级 结构 ,发 现 tRNA 缺 
Z DHU 臂 ,而 其 余 的 21 个 tRNA 基因 均 呈 现 出 典 
型 的 三 叶 草 结构 ,这 与 鳞 翅 目 中 目前 已 得 到 的 其 他 
昆虫 线粒体 基因 组 中 tRNA 基因 一 致 ( 王 菊 平 等 ， 
2012) ,也 与 很 多 后 生动 物 的 线粒体 基因 组 tRNA 基 
因 的 表现 ( 张 方 和 米 志 勇 , 1998) 相 一 致 。 在 小 红 珠 
绢 蝶 线 粒 体 基因 组 22 个 tRNA 基因 的 二 级 结构 中 








出 现 了 14 X (tRNA, tRNA®™, tRNA™, tRNA®, 
tRNA™, | tRNA — tRNA tRNA, 
tRNA™, tRNA™, iRNA?"99 tRNA™, (RNA T 
tRNA“) G = U 碱 基 错 配 现象 ,其 在 tRNA 结构 的 氮 
AENEA EE. J ROUTE, UK 5 e Al The 臂 均 
有 出 现 ; 此 外 ,出 现 1 对 (tRNA ) U = U 碱 基 错 
配 现象 ,其 分 布 于 反 密 码 子 臂 。 碱 基 错 配 现象 在 昆 
虫 中 较为 普遍 (Adams et al., 2000; Bae et al., 
2004.) , Yokobori 和 Paabo( 1995 ) 推测 造成 此 现象 的 
原因 可 能 是 由 于 昆虫 线粒体 基因 组 DNA 重组 功能 
的 欠缺 所 致 。 
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到 3 小 红 珠 绢 蝶 线 粒 体 基因 组 22 个 tRNA 基因 二 级 结构 
Fig. 3 Secondary structure for 22 tRNA genes from Parnassius nomion mtDNA 


2.5 


因 ,可 与 蛋白 质 


结合 


rRNA 基因 分 析 
rRNA 基因 是 一 类 相对 分 子 量 较 大 的 RNA 基 
形成 核糖 体 。 小 红 珠 绢 蝶 线 粒 
体 基因 组 中 包含 两 个 rRNA 基因 , 即 rrnL 和 rrnS, 分 
别 位 于 RNAL 和 tRNA 之 间 以 及 RNA“ #l 
A + 了 富 含 区 之 间 ( 图 1)。 两 个 基因 的 长 度 分 别 为 
1 348 bp (rmL) 和 797 bp (rms), AA H KES 
GenBank 中 鳞 翅 目 其 他 昆虫 相 比 , 都 属于 中 等 长 度 
(3X 3), rrnL 和 rrns 基因 的 A+T 含 量 分 别 为 


84. 4% 和 85.2% ,与 已 知 鳞 翅 目 昆虫 线粒体 基因 组 


( 表 5) 


o 


2.6 A«TEZIX 
A+T 语 含 区 为 线粒体 基因 组 中 的 非 编 码 序 列 ， 








中 rrnL 和 rrnS 基因 相似 , 均 含 有 较 高 的 A+T 含量 


即 为 不 能 转录 为 相应 的 mRNA ,从 而 不 能 指导 合成 
相应 蛋白 质 的 区 域 。 小 红 珠 绢 蝶 线 粒 体 基因 组 的 
A+T 富 含 区 位 于 rms 基因 与 tRNA 基因 之 间 ( 图 





1) ,其 长 度 为 483 bp, jp T X 28 WW PK RU 4 RE UR 
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之 间 。 小 红 珠 绢 蝶 线 粒 体 基因 组 A+T 寅 含 区 的 
A +T 含 量 为 93.0%( 表 4) ,高 于 其 线粒体 基因 组 其 
他 区 域 的 A + 了 TT 含量。 此 外 ,通过 与 其 他 鳞 翅 目 昆 
虫 比 对 分 析 发 现 ,小 红 珠 绢 蝶 线 粒 体 基因 组 的 A+T 
富 含 区 存在 以 下 结构 特点 :首先 ,rrns 基因 下 游 有 一 
段 由 “ATAGA"” 引 导 的 一 段 长 度 为 18 bp B SE T 
结构 (图 4) ,这 段 多 聚 了 结构 在 所 有 已 测定 的 鳞 翅 
目 昆虫 mtDNA 中 都 很 保守 ,长 度 从 12 bp( 宽 尾 风 蝶 
Agehana maraho ) 到 22 bp ( Hf 3& f Bombyx 
































rrnS--- 


mandarina) 不 等 (Kim et al., 2010) ,已 被 证 实 为 复 
制 起 点 。 其 次 ,在 小 红 珠 绢 蝶 的 A+T 富 含 区 序列 
中 存在 2 段 长 度 为 15 bp 的 重复 单元 ,重复 序列 之 
间 的 间隔 为 164 pp。 与 鳞 翅 目 其 他 昆虫 线粒体 基因 
组 A+ 了 T 寅 含 区 相 比较 , 小 红 珠 绢 蝶 该 区 存在 (AT)， 
结构 。 最 后 ,在 IRNA“ 的 上 游 可 见 一 个 17 bp 的 多 
聚 A 结构 ,被 认为 与 A+T 寅 含 区 的 转录 有 关 。 另 
外 ,小 红 珠 绢 蝶 线 粒 体 A +T 富 含 区 序列 中 散在 着 
一 些 长 短 不 一 的 串联 重复 单元 。 

















AAATTTTCTCATATAGATTTTTTTITITTTIITTIITATATTAATATATTTATTAA ATAATATAA 





复制 起 点 





The origin of replication 


TAAAATCTAAAATATTATTTTTTATACTAATATCAATATAAAAATTTGAATTTTTTTTTTTTT 
TTTTTTATTTTTAAAATTTTTTATAATAAATTTTTATTTTTTTTAATTTTTTAAAATAAAAAA 
AAAAAAAGTATATAAAATATTTAATATCAATTATTAAATATTGAAAAATTTCTTTTTTTTT 


15 bp 的 重复 单元 
15 bp repeat element 


TTTTCATAATTTGGAATAAAAACCCAAAATGGTTTTTTAAATTTTTTAAATTCAAAAAAA 
TTTTTTTAAATTAAAATATTTAAAAATTATATAAGTGTATTTAATAA ATTAAAATAAAATTT 
AAAAATATAAAATATTTAAAAAAAAAAAATATATATGCTACACCAAATTAAGGTATTTTT 


15 bp 的 重复 单元 











15 bp repeat element 


微 卫星 (ATn 单元 


Macrosatellite (AT)n element 


CTTTTATTTATTCCATTAACCCTTTTTTATTTTTTTTCTTATATAAAAAAAAAAAAAAAA A, 




















多 聚 A 结构 
Poly-A element 
-—IRNAP* 
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Fig. 4 Structure of A + T-rich region of Parnassius nomion mitochondrial genome 


2.7 蝶 类 高 级 阶 元 系统 发 育 关系 

以 美国 白 蛾 、 家 大 和 亚洲 玉米 螟 作为 外 群 ,基于 
28 种 蝴蝶 的 全 线粒体 基因 组 的 13 个 蛋白 质 编 码 基 
因 核 昔 酸 序列 建立 了 MP 树 和 贝 叶 斯 树 。 结 果 ( 图 
5 MIE 6) T78 , KR DL LP BUDE BEES JB ELA PR BJ 
RARR EPIRI DUREE BRE R Wl BE 
IRE, DPR SL BLUE El AI JK BI ËL BJ ZR 2ë 2 ET ORAS E 
近 , 而 与 凤 蝶 科 和 绢 蝶 类 物种 的 亲缘 关系 较 远 。 其 
中 眼 蝶 科 、 粉 蝶 科 和 灰 蝶 科 为 单 系 类 群 ,而 凤 蝶 科 内 
的 锯 凤 蝶 亚 科 首 先 与 绢 蝶 科 肾 为 一 文 后 青 与 风 蝶 亚 
科 的 物种 相 肉 。 结 果 表 明 , 绢 蝶 科 物种 与 锯 凤 蝶 亚 
科 物 种 的 亲缘 关系 最 近 , 这 一 结果 与 周 芜 先 生 将 绢 
蝶 类 物种 作为 独立 于 风 蝶 科 外 的 绢 蝶 科 处 理 和 锯 凤 
蝶 亚 科 与 凤 蝶 亚 科 亲缘 关系 更 近 的 观点 不 同 ;但 与 


























































































































武 春生 、 丁 亮 和 部 分 国外 专家 的 观点 一 致 ( Munroe , 
1961; Hancock, 1983; Igarashi, 1984; 武 春生 ， 
2001; 丁 亮 等 , 2007 ) ,该 结果 支持 将 绢 蝶 科 物种 归 
为 绢 蝶 亚 科 , 绢 蝶 亚 科 与 银 风 蝶 亚 科 及 凤 蝶 亚 科 统 
一 归 入 凤 蝶 科 内 , 且 绢 蝶 亚 科 与 句 凤 蝶 亚 科 的 亲缘 
关系 最 近 , 二 者 互 为 姐妹 群 。 由 于 测 线 粒 体 基 因 组 
全 序列 的 物种 有 限 ,目前 的 研究 仍 无 法 准确 确定 绢 
蝶 类 , 锯 凤 蝶 亚 科 和 凤 蝶 亚 科 的 分 类 地 位 及 其 之 间 
的 系统 发 育 关系 。 

在 风 蝶 科 绢 蝶 亚 科 内 ,小 红 珠 绢 蝶 与 其 他 5 种 
绢 蝶 科 物种 的 系统 发 育 关系 为 : (元 首 绢 蝶 P. 
cephalus + (君主 绢 蝶 P. imperator + ( fk |H 28 BÉ P. 
epaphus + 小 红 珠 绢 蝶 P. nomion) + ( 红 珠 绢 蝶 P. 
bremer + 阿波 罗 绢 蝶 P. apollo) ) ) 。 结 果 表 明 ,小 红 
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图 5 
Fig. 5 Bayesian phylogenetic tree of butterflies based on the sequences of 13 protein-coding genes of mitochondrial genome 
物种 中 文 名 和 GenBank 登录 号 见 表 2。 节 点 处 的 数值 为 该 节点 的 后 验 概 率 。 图 


Bayesian posterior probabilities are indicated at each node. The same for Fig. 6. 


Bombyx mori 
Hyphantria cunea 
Ostrinia furnacalis 
Parnara guttata 
Heteropterus morpheus 
Libythea celtis 
Danaus chrysippus 
Callerebia suroia 
Lethe dura 
Melanargia asiatica 
Acraea issoria 
Apatura ilia 
Abisara fylloides 
Lycaena phlaeas 
Japonica lutea 
— Cupido argiades 
Leptidea morsei 
Eurema hecabe 
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Pieris melete 
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Papilio protenor 
Papilio glaucus 
Papilio maraho 
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Parnassius imperator 
Parnassius bremeri 
Parnassius apollo 
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基于 线粒体 基因 组 13 个 蛋白 编码 基 





因 序列 构建 的 蝶 类 贝 叶 

















6 同 。Species name and GenBank accession number see Table 2. 





珠 绢 蝶 与 依 帕 绢 蝶 P. epaphus ZIRA HR P. bremer 
和 阿波 罗 绢 蝶 P. apollo 的 亲缘 关系 较 近 ,其 中 小 红 
珠 绢 蝶 P. nomion 与 依 帕 绢 蝶 P. epaphus 的 亲缘 关 
系 最 近 。 





3 讨论 与 结论 











组 全 序列 具有 很 高 的 A + T 偏向 性 ,验证 了 昆虫 线 
粒 体 基因 组 中 存在 A + 了 含量 偏向 性 的 特征 。 此 
外 ,通过 计算 小 红 珠 绢 蝶 线 粒 体 基因 组 全 序列 的 AT 
Ai ERI GC 偏 斜 ,发 现 其 序列 中 有 较 多 的 T 和 C, 这 
与 风 蝶 总 科 中 的 大 卫 绢 峡 蝶 C，buddha， 蔚 麻 珍 蝶 
Acraea issoria, X] JK WE Coreana raphaelis， 仁 眼 蝶 
Hipparchia autonoe， 红 珠 绢 蝶 Parnassius bremeri 和 



























































小 红 珠 绢 蝶 线 粒 体 基 因 组 中 ,37 个 基因 的 排列 KREKAR Teinopalpus aureus 的 相 一 致 ( 夏 靖 等 ， 
顺序 和 编码 的 方向 与 已 报道 的 大 紫 峡 蝶 和 大 卫 绢 峡 — 2011). 











WE C. buddha 相 一 致 (Kim et al., 2006; Hong et al., 
2009; Kim MI et al., 2009; Hu et al., 2010; Kim et 
al., 2010) , 且 未 发 现 基 因 重 排 和 缺失 现象 。 通 过 线 
粒 体 基 因 组 核 音 酸 组 成 分 析 发 现 小 红 珠 绢 蝶 线 粒 体 
基因 组 13 个 蛋白 质 编 码 基 因 、2 个 rRNA 基因 和 
A+T 富 含 区 中 A + 了 碱 基数 占 总 碱 基数 的 比例 明显 
高 于 G+C 的 比例 ,表明 小 红 珠 绢 蝶 的 线粒体 基因 









































与 其 他 昆虫 相 比 , 鳞 翅 目 昆 虫 的 线粒体 基因 组 
的 A+T 寅 含 区 相对 保守 。 鳝 埃 目 昆虫 A+T 富 含 
区 中 存在 (AT), 的 结构 ,小 红 珠 绢 蝶 为 (AT): ,大 紫 
WEE 7g (AT), KEARI Jg (AT), EJ JK HE ON 
CAT); , IIRA RO CAT) o AZ N CAT) u (Kim 
et al., 2006; Hong et al., 2009 ; Kim MI et al., 2009; 
夏 靖 等 , 2011; 王 菊 平 等 , 2012) 。 另 外 , 小 红 珠 绢 
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Ostrinia furnacalis 
53 Bombyx mori 
L= aut ais 
100 Pamara guttata 
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图 6 基于 线粒体 基因 组 13 AE AA 36442 RF PP EE B WI AA o K i 24] 22 2: EET 








Fig. 6 Maximum parsimonious phylogenetic tree of butterflies based on the sequences of 


13 protein-coding genes of mitochondrial genome 


分 支 节 点 的 数值 表示 1 000 次 内 部 分 支 检 验 置信 值 。Bootstrap values for 1 000 replicates are indicated at each node. 








蝶 线 粒 体 A + 了 T 富 含 区 序列 中 散在 着 一 些 长 短 不 
的 串联 重复 单元 ,但 在 与 其 近 缘 物种 非 编码 区 的 比 
较 中 并 未 发 现 相似 的 保守 序列 区 ,因此 ,认为 此 At+ 
T 富 含 区 中 多 数 序列 对 线粒体 基因 组 无 调控 功能 
( Broughton et al., 2001) 。 

本 研究 测定 了 小 红 珠 绢 蝶 线 粒 体 全 基因 组 序列 
并 构建 了 包含 弄 蝶 总 科 和 凤 蝶 总 科 的 28 种 蝴蝶 的 
系统 发 育 树 ,探究 蝶 类 高 级 阶 元 的 系统 发 育 关 系 , 分 
析 结 果 与 武 春生 、 丁 亮 及 部 分 国外 专家 的 观点 一 致 
(Munroe, 1961; Hancock, 1983; Igarashi, 1984; 武 
春生 , 2001; p 2626, 2007) 。 该 结果 支持 将 绢 蝶 科 
物种 归 为 绢 蝶 亚 科 , 绢 蝶 亚 科 与 锯 凤 蝶 亚 科 及 凤 蝶 
亚 科 统 一 归 入 凤 蝶 科 内 , 且 绢 蝶 亚 科 与 饮 凤 蝶 亚 科 
的 亲缘 关系 最 近 , 二 者 互 为 姐妹 群 。 在 凤 蝶 科 绢 蝶 
亚 科 内 ,小 红 珠 绢 蝶 P. nomion 与 依 帕 绢 蝶 P. 
epaphus 的 亲缘 关系 最 近 。 这 一 结果 与 传统 的 形态 

































































学 分 类 结果 基本 一 致 。 本 文 结果 不 仅 为 明确 绢 蝶 科 
和 凤 蝶 科 的 分 类 地 位 提供 了 重要 的 分 子 依 据 , 也 为 
后 续 进 一 步 对 绢 蝶 科 物种 进行 线粒体 基因 组 测定 和 
分 析 , 澄 清 绢 蝶 科 各 物种 之 间 的 系统 发 育 关 系 葛 害 
了 分 子 基 础 。 
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